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ABSTRACT 

In this paper we present a theoretical and numerical study of the growth of 
linear disturbances in the high-Reynolds-number laminar compressible wake behind a 
flat plate which is aligned with a uniform stream. No ad hoc assumptions are made 
as to the nature of the undisturbed flow (in contrast to previous investigations) but 
instead the theory is developed rationally by use of proper wake-profiles which satisfy 
the steady equations of motion. The initial growth of near wake perturbations is 
governed by the compressible Rayleigh equation which is studied analytically for 
long- and short-waves. These solutions emphasize the asymptotic structures involved 
and provide a rational basis for a nonlinear development. The evolution of arbitrary 
wavelength perturbations is addressed numerically and spatial stability solutions are 
presented that account for the relative importance of the different physical mechan- 
isms present, such as three-dimensionality, increasing Mach numbers and the non- 
parallel nature of the mean flow. Our findings indicate that for low enough (sub- 
sonic) Mach numbers, there exists a region of absolute instability very close to the 
trailing-edge with the majority of the wake being convectively unstable. At higher 
Mach numbers (but still not large - hypersonic) the absolute instability region seems 
to disappear and the maximum available growth-rates decrease considerably. Three- 
dimensional perturbations provide the highest spatial growth-rates. 


1 This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NAS1-18605 while the author was in residence at the Institute for Computer Applications in 
Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665. 


- 2 - 


1. Introduction. 

In recent years there has been an increased interest in compressible shear layer flows due to the 
need in development of hypersonic propulsion systems for trans-atmospheric flight. The major problem 
that faces these projects, is the degree of stabilization that compressible mixing layers exhibit as the 
Mach number increases. This in turn means that the high degree of mixing that is desirable before igni- 
tion (in the scram-jet engine for example) is not available inside the system. A quantitative understand- 
ing of the instability mechanisms present is of much value, therefore. Some techniques have been sug- 
gested to alleviate this problem, by Kumar et al. (1987). 

Most of the effort in understanding the stability of shear layers has been centered around the fully 
developed monotonic compressible shear layer which is typically modeled by a hyperbolic tangent dis- 
tribution, together with the corresponding temperature profile (see for example Gropengiesser (1969), 
Jackson and Grosch (1989a,b), Ragab and Wu (1988)). In practical applications, such as those sited 
above, the mixing occurs when two streams come together at the trailing edge of a splitter plate to 
form a wake. If the free-stream velocities of the two components are different (usually this is the case, 
but in many cases they are numerically close) the tanh shear layer forms the far wake development of 
the steady equations of motion. There is a region closer to the trailing edge, where the unperturbed 
flow has a definite wake component. It is the purpose of this paper to analyze the instability mechan- 
isms of this wake region. In order to fix matters we will consider the case of a plane wake where the 
free-stream values of the two oncoming fluid components are equal. 

Incompressible wakes have been studied extensively for the existence of unstable waves through 
numerical and analytical solutions of the classical Rayleigh stability equation. Most early studies 
assume an ad hoc basic flow in the form of a gaussian or sech 2 velocity distribution (Hollingdale 
(1940); McKoen (1957); Sato and Kuriki (1961); Mattingly and Criminale (1972)). It is clear that such 
basic flow does not satisfy the equations of motion and at best provides the correct characteristics of 
the far-wake. Further more, the flow characteristics in the near wake are fairly complicated and involve 
the double structure described in the seminal work of Goldstein (1930), where there is a displaced 
outer boundary-layer region together with a thinner region including the wake center-line through 
which the abrupt change of boundary conditions (from no-slip upstream, to symmetry downstream) is 
accommodated. Papngeorgiou and Smith (1989) (referred to as PS1) have considered incompressible 
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wake-stability, both analytically and numerically, using a correct basic flow that satisfies the equations 
of motion for large Reynolds numbers (the wake-boundary-layer equations of Goldstein). The non- 
linear development of the instability is analyzed in an accompanying paper, Papageorgiou and Smith 
(1988) (referred to as PS2). 

A very interesting and physically important aspect of shear layer flows is the possibility of con- 
vective or absolute instabilities. For an account of the ideas involved in shear layer stability see Huerre 
and Monkewitz (1985), Huerre (1987). This aspect of instability has been studied by Betchov and 
Criminale (1966) and more recently by Hultgren and Aggarwal (1987), for a two-parameter family of 
model basic states. More specifically the basic flow is assumed to have the form u = 1 - a e -by , where 
0 < a < 1 , b > 0 or an equivalent scch 2 distribution. Hultgren and Aggarwal for example, find that 
the flow becomes absolutely unstable when a £ 0.943. Such a profile clearly corresponds to a station 
very close to the trailing edge as we briefly demonstrate. The center-line velocity of the proposed 
absolute-instability profiles is at most 0.057 in non-dimensional terms. Using Goldstein’s results (see 
also Smith (1984)) we know that the velocity in the near-wake at the center-line is given by 
1 .6 1 1 x i/3 A, 2/3 where X is the Blasius skin-friction and is equal to 0.33206... . Equating the value 0.057 
with this, gives the value for the non-dimensional x-station to be x = 0.0004017, or £ = x ,/3 = 0.0738. 
The basic flow at this point does not seem to be representable by the gaussian model (see our results 
section) and so some serious objections may be raised as to the rational prediction of 
convective/absolute change-overs by use of such models. The trends seem to be correct, however, in 
comparison to the results of PS1 which suggest an absolute instability region very near the trailing 
edge. The important question that arises in supersonic wakes is, besides the possibility of absolute ins- 
tability regions per se, the characteristics of an absolute instability region as the Mach number 
increases. This is addressed numerically for two Mach numbers that are representative of the different 
physics that is involved; these are 0.1 and 3.0 corresponding to compressible subsonic and supersonic 
flow respectively. It is found that the subsonic flow contains a region of absolute instability whereas 
the supersonic one seems not to. The streamwise station at which the change-over from convective to 
absolute instability takes place for M=0.1 is less than 0.000001 in non-dimensional terms. 

The purpose of the present work is to extend the incompressible study of PS1 to compressible 
flows of arbitrary Mach numbers. Our approach is that of classical linear stability theory that covers 
both two- and three-dimensional disturbances. In Section 2 we briefly describe the basic flow and 
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formulate the linear stability problem that leads to the compressible Rayleigh equation. Section 3 is 
used to analyze the stability of short-wave disturbances in the very near wake where the compressible 
analogue of the Goldstein double structure provides an asymptotic description for the unperturbed flow. 
In Section 4 we consider the behavior of long-waves (these waves are long compared to the boundary- 
layer thickness but not as long as the plate length-see Section 2 also). Section 5 discusses numerical 
solutions of the Rayleigh equation using the nominally exact unperturbed profiles that are in turn com- 
puted numerically. In Section 6 we present numerical results for a range of streamwise stations, Mach 
numbers and angles of wave propagation, and finally we draw some conclusions. 

Throughout the paper we use a cartesian frame of coordinates which is fixed at the trailing edge 
x = 0. ’I he flow quantities p,u,v,w,p,T,(.t denote non-dimensional density, velocities in the x, 
y and z directions, pressure, temperature lihd viscosity coefficient respectively, where non- 
dimensionalizations are made with respect to free-stream values which are denoted by a subscript °o. 

U„L 

The Reynolds number is R c = — — where L is the plate length, the Mach number is M = U^a where 

a is the sound-speed and the Prandtl number is Pr = p„ Cp/k where c p is the specific heat at constant 
pressure and k is the coefficient of thermal diffusivity. R c is taken to be asymptotically large 
throughout, while is arbitrary within the context of validity of the equations of motion. The 
governing equations used are the 2D steady and pressure-free (classical) compressible boundary-layer 
equations for the calculation of the basic flow (bars are used to distinguish the basic flow), or the 
unsteady compressible 3D Euler equations for the evolution of wavy disturbances. We give these here 
for future reference. 

J) ( u u„ + V TTy ) = ( JT u y ) y , (1.1a) 

jT(uT x + 7T y ) = -L ( |I T y ) y + (y- I )M 2 frir y 2 ", (1.1b) 

( ptr ) x + 1 P v ) y = o , (i.ic) 

PT=1 . (l.ld) 

The boundary conditions are 


u y ( X , 0 ) = Ty ( X , 0 ) = 0 , u(x,oo) = T(x,oo) =r 1 . (1.1 BC) 


The 3D Euler equations are used in the form 
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p ( U t + U U x + V U y + W u z ) = - Px , 

, 1 

P ( v, + U v x + V v y + w v z ) = - Py , 

. 1 

p(w, + uw, + vw y +ww I ) = — Px , 

yM 

p ( T, + u T, + v T y + w T, ) = ( p, + u p x + V p y + w p, ) 

P, + ( p U ) X + ( P V )y + ( P W )„ . = 0 , 

P = P T (Y - 1) M 2 . 


(1.2a) 

(12b) 

(1.2c) 

( 1 .2d) 

(1.2e) 

(1.20 


2. Basic flow and formulation of the linear stability problem. 

The formulation that leads to the compressible Rayleigh equation essentially follows that of the 
incompressible analogue, as described in PS1 for example. It is instructive, however, to bear in mind 
the scales for which the Rayleigh equation is valid in the case of a planar wake. The arguments of 
PS1 carry over to the present problem also (the important physical mechanism is the asymptotically 
large Reynolds number which is common to both flows), and the instability develops on a local region 
of horizontal and lateral extent of order R" 1 ' 2 ; the time scale of the instability is also of order R c " 1/2 . 
These scalings insure the inviscid development of wavy disturbances of the viscous underlying basic 
flow. Further more, as is usual in Rayleigh-type evolutions, the scales are such that the x-variation of 
the basic flow as well as its lateral component, are higher order effects. In the present problem the 
streamwise component of the basic flow varies on an order one x-scalc while the lateral component is 
of order R' ,/2 . Viscosity becomes important when the wavelength of the perturbation becomes of 
order unity as compared to the wake-thickness. For more details of the scaling analysis see PS1 and 
PS2. 

Before the introduction of a perturbation it is useful to summarize the characteristics of the under- 
lying unperturbed state. This is described in detail by Papageorgiou (1989) and we sketch the results 
here. The governing equations are (1.1 a-d) together with the boundary conditions (1.1BC). Next we 
assume a model fluid and a linear viscosity law, i.e. the Prandtl number Pr = 1 and p = C T 
(Chapman’s law, C is a constant that is scaled out of the problem and the equation p = T is used ). 
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y. 

The Dorodnitsyn-IIowarth transformation Y = is employed, from which it is seen that Crocco’s 

equation (see Stewartson (1964)) becomes 

d 2 T . . . .,2 

— = -( y - 1 ) M z . 

fTu 

It follows that T involves two constants of integration which can be found by matching with the 
upstream flow. Two physically distinct cases are given, (i) constant wall temperature, T w say, on the 
plate, and (ii) insulated walls (i.e. the temperature gradient is zero on the plate upstream). In case (i) 
the result is 

T(x, Y)= 1 + -|- (7-1 )M 2 (u — u a ) + (T w - 1 )( 1 — u) , (2.1) 

while in case (ii) we have 

T( x , Y ) = 1 + -y ( y-1 ) M 2 ( 1 - u 2 ) . (2.2) 

In order to solve for the velocity u(x,Y), conditions must also be given at the outset x = 0 since the 
problem is parabolic in x. The nonlinear system that gives the undisturbed flow is 

u = tj/y , 11 U x — \j/ x Uy = Uyy , (2.3) 

holding for x > 0, together with the boundary conditions 

u( x , oo ) = l , \j7( x , 0 ) = u Y ( x , 0 ) = 0 , y( 0+ , Y ) = \j/ B ( Y ) , (2.3BC) 

where Yi/ represents the Blasius velocity profile given by 

Vb'" + \ Vb Vb" = 0 , Yd(0> = Vd'(0) = Vb'(~>) -1=0; Yb"(0) = ; X = 0.332... 

We will be mainly concerned with adiabatic wall conditions but the results can be extended to constant 
wall temperatures also. 

With the basic flow assumed known, a perturbation of size Ei is introduced so that the total flow 
has the form 

( u , v , w , p , T , p ) = ( u , 0 , 0 , jT , f , 0 ) + f , ( u , <j> , w , p , t , p ) e‘ ( ra x + p 7 ~ ‘ > , 
where the tilde functions depend on y alone. Substitution into (1.2a-f) (we are anticipating the inviscid 
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result here; equivalently a substitution into the Navier-Stokes equations together with the instability 
scales described earlier yields the same result), and linearization with respect to e provides a linear sys- 
tem of five ordinary differential equations. We will work with the transformed variable Y rather than 
the physical one and will use the following forms of the compressible Rayleigh equation obtained by 
elimination of four unknowns and dropping of the tildes : 

p" _ J^L p'-T 2 [ o 2 + [) 2 - ( tr-c ) 2 M 2 | p = 0 , (2.4) 

u - c T 


d ( u - c ) 4>' - XT' 4> 

dY G T 


- a 2 ( u - c ) <J> = 0 , 


(2.5) 


where primes denote Y-derivatives and G is given by 

G = — + P - T - M 2 ( u - c f . 
a 2 

Equations (2.4) and (2.5) describe the evolution of perturbation pressure and normal velocity respec- 
tively; (2.5) was used by Gropengiesser (1969) in his study of shear layers. There is a choice in boun- 
dary conditions at the wake center-line Y = 0, due to the symmetry of the basic flow. We choose to 
study two distinct wake-modes, anti-symmetric ones (mode 1) and symmetric ones (mode II). The con- 
vention used here is that of incompressible wake stability whereby the perturbation is termed anti- 
symmetric/symmetric if the streamwise velocity perturbation, u, is anti-symmetric/symmetric with 
respect to Y=0. It is generally found that the anti-symmetric ones are the most dangerous ones. Experi- 
mental evidence lends weight to this finding through the well-known Karman vortex-street found in 
wakes. The boundary conditions used are, therefore, 

Mode I p(Y=0) = 0 , p(Y=*») = 0 , Mode II p'(Y=0) = 0 , p(Y=») = 0 , (2.4BC) 

Mode I f(Y=0) = 0 , $(Y=°°) = 0 , Mode II <t>(Y=0) = 0 , <J>(Y=«) = 0 . (2.5BC) 

Here we have assumed that the disturbance decays at infinity. This assumption is reasonable since, in 
contrast to shear layers, there does not exist a supersonic region that can support neutrally stable oscil- 
lafory outgoing waves (see comments of Section 6 also; for a discussion of supersonic shear waves see 
Jackson and Grosch (1989a)). The eigenvalue problems defined by (2.4), (2.4BC) or (2.5), (2.5BC) 
are equivalent and must, in general, be addressed numerically. This is done in Section 5, but before 
discussing numerical solutions we consider two important limits which are amenable to analysis. These 
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limits are short- and long-waves which in 2D temporal theory are characterized by a — > <*> and a — » 0 
respectively (see PS1). It is worth noting that short- and long-wave analyses only make sense for tem- 
poral theories, but as it turns out the growth rates are asymptotically small in both cases thereby mak- 
ing the results readily applicable to spatial stability via the Gaster (1968) group velocity transformation. 
The analysis of the following sections considers two-dimensional waves and these results are extendi- 
ble to three-dimensional ones after minor modifications. 

3. Short waves, a — » oo, 

In this section we construct the asymptotic behavior of the solutions in the limit of large a. The 
physical structure of these waves is that they have length-scales which are small compared to the wake 
thickness. Further more we concentrate on the near-wake where the basic flow is described by the 
Goldstein double structure (see Papageorgiou (1989) for example). Our analysis can be viewed in 
terms of the double limit x > 0 and ot — > °°, The ideas and structures involved are similar to the 
incompressible case studied in PS1, but the growth rates we compute here are larger than their 
incompressible counterparts. 

The starting point is the Rayleigh equation (2.5) which is considered for 2D anti -symmetric 
(mode I) disturbances - equivalently the pressure equation can be used, as well as mode II waves but 
as the subsequent numerical results indicate, mode I waves are more unstable (see comments in the 
introduction also). Inspection suggests the following expansions in the re-scaled region Y = a -1 rj with 
r \ = 0 ( 1 ) : 

4 > = < t > 0 + a " 1 < J >, + a ' 2 <(> 2 + • • • , (3.1a) 

u = ± a ' 1 X i\ + a " 4 X t r | 4 + , (3.1b) 

c = c 0 + a -1 c, + of* c 2 + of 3 c 3 + • • • . (3.1c) 

For adiabatic walls the temperature expands as 

T = T 0 - cf 2 Ti 2 T 2 + • ■ • , T 0 = 1 + ^(y-l)M 2 , T 2 = y(y-l)\ 2 M 2 (3. Id) 

The constant X in (3.1b) is the skin-friction of the Dlasius flow discussed in the previous section, and 
= -(1/48) X 2 . The expansion (3.1b) reflects the fact that we are outside the Goldstein layer where 
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the basic flow is essentially the displaced Blasius flow. In the analysis that follows the Goldstein layer 
(of thickness x I/3 ) is thin compared to the wavelength 2tt/a; the ordering a -1 :*> x 1/3 is therefore 
observed. Physically this corresponds to the development of waves which are short compared to the 
boundary-layer thickness but not as short as the Goldstein layer thickness. It can also be seen that 
given any large a , the position x can be chosen to be small enough so as to validate the above order- 
ing (this is achievable as long as x » R,T 3/8 which means that we are still outside the triple-deck). 

In view of the above, the Goldstein layer can be regarded as a sharp interface across which the 
perturbation pressure and normal velocity must be continuous. This means that the following conditions 
hold to all orders in a -1 : 

[<H = 0 , [(u-c)<t> n -u<|>] = 0 , (3.2a, b) 

where the square bracket denotes the jump across r| = 0. In view of these physical constraints, we find 

that .for a non-trivial solution the leading order phase-speed vanishes, i.e. 

i 

I 

; c 0 = o . 

is not surprising in view of the fact that the basic flow has a maximum streamwise component of 
a -1 with which it can convect disturbances downstream. Substitution of (3.1a-d) into (2.5) gives 
■flowing solution to leading order 

4>o = A 0 e~ + T ° n . (3.3) 


± refers to the solution for fl > 0 and t| < 0 respectively. Condition (3.2a) has been used to 
le same constant A 0 in the upper and lower half planes respectively. When the leading order 
n (3.3) is used in the pressure continuity condition (3.2b), and 0(a _I ) are balanced, we obtain 
ding contribution to the phase-speed to be 



(3.4) 


the compressible analogue of the result in PS1. It can be seen from this expression that the 
leed of the perturbation decreases with increasing Mach number (cf. (3. Id)) and in fact decays 
call^y in the hypersonic limit. Next we find 4> 2 = 0 , c 2 = 0 (see below). This implies that the 
ives constructed here are almost neutral and have growth-rates of 0(cT 3 ) at most. This is in 
to the incompressible problem which yields O(oc~ 4 ) growth-rates (cf. PS1). The complex 
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coefficient c 3 is computed by construction of uniformly valid asymptotic expansions which are matched 
across the critical layers and the central Goldstein layer to satisfy the boundary conditions and so pro- 
vide the correction to the eigenvalue problem. 

Before proceeding to the higher order approximations, it is useful to introduce a new vertical 
coordinate given by ^ = T 0 T|; this simplifies the algebraic expressions since the critical layers now 
occur at £ = ± 1 (see below). The next corrections to the eigenfunction satisfy the following equa- 
tions, then 


( ± £ - 1 ) ( <J>i^ - <h ) = 0 , 4>,(±«>) = 0 , 

<t>2f£ - 4>i = ± 2 S (y $ -+1) e-* * - S(-£ 2 -+ 7% + 1) e^ * 


(±£- 1) 


^ , <j) 2 (± 00 ) = 0 , 


where as Wore the upper sign is to be taken in the region £ > 0 and the lower one in £ < 0. The con- 
stant S appearing in 3.6 is given by 

„ * A 0 M 2 


The solution for <J>, is identical to that of <t> 0 and so, without loss of generality, it is taken to be zero by 
re-normalization of A 0 . As mentioned earlier, it can be seen that (3.6) has singularities at the levels 
^ = h these singularities are critical layers and arise because the phase-speed is neutral to leading 
order (cf. (3.4) above). The solution for (j >2 then splits up into four distinct regions (cf. Figure 1 in 
PS1). The four regions 1 < £ < , 0 < £ < 1 , -1 < £ < 0 , -~ < ^ < -1 are denoted by 

(1), (2), (3) and (4) respectively and the following solutions are found 

♦1° * e* [ B, - S (“^ 3 + j (Y+l) 4 2 + (Y~|) ^ + (y-2) in(jj - 1) 


+ S e^J [ -±tf + (Y f l) + (y- 2) + 1 e^' d§, , (3.7a) 

z 1 Sl~» 

»f> = e* [ Bj - S (-jti 1 + y (nl) Z 1 + (Y-y) Z + (y-2) In|5 - i| 

+ A, e 5 + S c E |( -yl^, 2 + (yt-l) 5, + (y-y) + le"’ 5 ' d?, . (3.7b) 
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<i>i 3) = e* [ B 3 + S (--^ 3 - T (YH) ¥ + OMr* ^ “ ( ^ 2) ,n( ^ + ]) 

6 2 2 

+ A, e« - S e-5 f [ -{rt? - (Yf I) 5, + <Y-f ) - ^ I e*‘ «■ • OK) 

$l a) = e 5 [ B„ + S (-{rt 3 - 4- (fU) t? + (r4> 5 - (Y-2) !"I5 + 1 1 

6 2 2 

- s e-t J I -{rtf - mi) 5, + (y-b - I «*' 15, • < 3 - 7d ) 

In matching solutions across the critical levels we assume that the constants in the various regions are 
suitably connected so as to render our inviscid solution as the zero viscosity limit of the Orr- 
Sommerfeld equation. This amounts to picking up the correct branch from the complex logarithmic 
dependence as has been explained by Lin (1955), Lees and Lin (1946) (see also Drazin and Reid 
(1981)). The net effect, therefore, is a phase-jump of -\k sgn(if(£ = ±1) ) as a critical level is crossed 
from above (the temperature basic state is passive in this because it is positive throughout the domain). 
Another set of connection equations is obtained from the continuity of <J> across the levels £ = ±1. The 
analysis is fairly lengthy so we just give the results here. The following four algebraic equations are 


obtained that connect the six unknowns Aj , A 2 , B, , B 2 , B 3 , B 4 : 

B, -B 2 -e 2 A 2 = e(p 3 -p,) , (3.8a) 

B, - B 2 +e 2 A 2 = p 4 - 27tiS(y-2) , (3.8b) 

B 3 - B 4 4e 2 A 3 = e (p 7 - p 5 ) , (3.8c) 

B 3 - B 4 -e 2 A 3 = -e (p 6 + p 8 ) + 2/tiS(y-2) , (3.8d) 


where Pi , p 3 , p 4 , p 3 , P6 > P7 » Ps are rea * constants that arise from the limiting values of individual 
solutions as the critical layers are approached from above or below. The exact effect of these constants 
is of secondary importance to our purpose here since they are real valued and do not contribute 
towards the instability. The complex terms (phase-jumps from the critical layers) are explicitly included 
in the equations (3.8a-d) above. 

The remaining two algebraic equations that enable complete solution for <() 2 are obtained by con- 
sideration of the Goldstein layer dynamics. We point out again that the scales involved are essentially 
those of the incompressible problem of PS1; the compressibility effects, however, introduce 
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perturbation dynamics at a lower order. The Goldstein layer has thickness e = x I/3 « 1. A rescaling of 
the Dorodnotsyn-Howarth variable is required, therefore, and the order one variable Z is introduced, 
given by Y = e Z. Consideration of the outer solutions for small rj (or £ ) suggests that the eigenfunc- 
tion inside the Goldstein layer proceeds in powers of ea (this comes from the expansion of <{> 0 ). The 
next outer correction is a -2 <t> 2 and as we have shown above <j> 2 contains an order-one displacement 
Further more it is clear that instability will arise from the contribution of <j) 2 (it contains the complex 
phase-jumps) and so we chose to structure our expansion so that the effects of <}> 2 come in at second 
order. This requires the balance ea = a -2 or e = a -3 . In view of the definition of e this scaling 
implies an asymptotic wavelength of order x as x tends to zero. The validity of the solutions can 
be improved by a double expansion in powers of a -1 and ea, whereby the asymptotic condition is 
ea <c: 1. Both methods give the same final result (as was found for the incompressible problem of PS1 
also) and for brevity we present the former analysis. 

Using a -1 as the perturbation parameter the Goldstein layer is given by Y = a -3 Z and the 


expansions are found to be \ 

♦o(Z) = d> 0 (Z) + a _2 0,(Z) + a“ 4 <I> 2 (Z) + • • • , (3.9a) 

u = a 3 F(Z) + 0(a~ 12 ) F(Z) - ± X ( |Z| + A G ) as Z->i« , (3.9b) 

T = T o - y cT 6 (y-l)M 2 F 2 + 0(a-15) , (3.9c) 

c = a~' c, + a“ 3 c 3 + • • • , (3.9d) 


where the function F(Z) is the self-similar near-wake solution found by Goldstein (1930); the constant 
A g represents the positive displacement of the upstream boundary-layers. The three leading order 
solutions are then found to be i'-r 

<I> 0 = A 0 , (3.10a) 

<I’l = F(Z) + A, Z + A 2 , (3.10b) 

c t 

<h' 2 = ( ) F'(Z) + T 0 2 A 0 Z - ( Z F'(Z) - F(Z) ) . (3. 10c) 

cf c i c, 


where A, , A 2 are constants to be determined from matching. Note that the solutions above reflect the 
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fact that we are considering mode I waves. With (3.10a-c) at hand, it is fairly straight-forward to match 
with the outer solutions as |Z| -> « and |^| -* 0 respectively. This implies that Aj = 0 and yields 
the two required algebraic equations: 


A2 + B 2 — A3 + B3 
Aj + A3 — — B3 = 


, t a 0 a.a g ^ 

( = A 2 ) ♦ 

C 1 

(3.11a) 

2k A 0 C ) A 2 . 

T 0 1 c, 2 c, • 

(3.11b) 


Solution for the six unknown constants yields the second order correction to the wave-speed c to be 


c 3 


k A _ 4?tiS(y-2) 
2T 0 A 0 [ e 2 


(3.12) 


In (3.12) above, A is a real constant which we need not give here. The imaginary part of c 3 is given 
by (using the definition of S) 

Cr = 27i(Y- 2 )X 3 M 2 (3.13) 

e 2 [1 + ( l/2)(y— 1)M 2 ] 4 

Since y < 2 (in fact y = 1.4 ) the imaginary part of c 3 is positive and so the flow is linearly unstable. 
This last remark follows from the correct identification of physically meaningful zero viscosity limit 
solutions of the Orr-Sommerfeld equation (cf. Lin (1955)). Another point worth noting, is that the 
growth rate in the hypersonic limit behaves like as M — > as opposed to the M decay of the 
leading order phase-speed C|. This means that the correction oT 3 c 3 always remains asymptotically 
smaller than of'cj as the Mach number increases without bound for any functional dependence of 
Mach number on a. 


4. Long waves, a — * 0. 

The intuitive approach to long-wavelength perturbations indicates that the exact details of the 
basic flow near the center-line are not that important, in obtaining leading order results at least, and the 
solutions depend on the basic flow characteristics at infinity. This was pointed out and utilized by 
Drazin and Howard (1966) in the description of several "broken-line" profiles, the results of which 
extend to wakes also. Originally the analysis described incompressible flows, but was later extended to 
include compressible ones also (Gill and Drazin (1965), Blumen et al. (1975)). The approach 
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essentially depends on obtaining solutions in the upper and lower half-planes respectively, and using 
their Wronskian to fix the eigenvalues. In the present work we present a different approach which is 
based on the scalings and structures of the various asymptotic regions involved. This is useful in that 
it motivates scalings for low-frequency/long-wavelength nonlinear stability also. The incompressible 
problem is described in PS1, and its low-frequency nonlinear stability analysis is given in PS2. The 
asymptotic limit a — » 0 is therefore considered and as in Section 3 we choose to work with equation 
(2.5) for the perturbation normal velocity. Physically, we arc considering instability scales that arc 
long compared to the wake thickness but still short compared to the plate length. Further more mode I 
disturbances are considered since they generally provide higher growth rates (see earlier comments). 

A region Y = 0(1) is considered first (region I), where the appropriate expansions are found to 

be : 


<|) = $0 + a'% + a<j> 2 + + a 2 <J> 4 + • ■ ■ , (4. la) 

1 n 

c = c 0 + a' % + a.c 2 + • ■ • . (4.1b) 

The function G that appears in (2.5) also expands as follows : 

G = G 0 + a 1/2 Gj + ■ • • , 

where 


G 0 = T(Y) - M 2 (u-c 0 ) 2 . (4.1c) 

Let us now define the following linear operator L : 

L(.) = (u-Co)-^- - vr ( ) • (4.2) 

Substitution of (4. la-c) into (2.5) gives the following four leading order problems : 



77- (LOhH^'o) 
Go 


77" (L(<|> 2 ) _c l < l> , l~< : 2 ( i >, o) 
Go 


* 

(4.3a) 

= 0 , 

(4.3b) 

1 = 0 , 

(4.3c) 
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d 

dY 


77 " (L(<t>3) _c l4 >, 2 -c 2 < t ), l _c 3 < t >, o) 
Go 


d 

dY 


77- (L(d>4)— c id >, 3 — c 2 < t ), 2 — c 3 < ? >, l -_c 4 < ( >, o) 

Go 


= 0 , 

(4.3d) 

(u-Co) A 
J - <Po ■ 

(4.3e) 


The appropriate solution of (4.3a) that is both even and bounded for large Y is the displacement solu- 
tion 

to = Ao(u-c 0 ) • (4 ' 4) 

and to ensure decay to zero at infinity the leading order wave-speed is 

c 0 = 1 . 

It is not surprising to find, therefore, that long-waves are almost neutral and have phase-speeds approx- 
imately equal to that of the free-stream. Next we find 

$,=- 0^0 , <t> 2 = - c 2 A 0 , <t>3 = “ c 3 A 0 • 

Substitution of (4.4) into (4.3e) gives the following solution for <j> 4 : 


Y G 0 ( y] )dy, 


„ t G 0 y, dy, ( (u-1 z dy 2 , J „ f 

= A 0 (u-l) I 7— I = +d 4 (u-l)J— 2 

i (u-i) 2 L t b (u-i) 2 


- C 4 A 0 , 


(4.5) 


where d 4 is given by 

d4 = A 0 p^dY , (4.6) 

since mode I disturbances are being analyzed and so <t>prim 4 (0) = 0. It is the boundary condition (4.6) 
that provides the solution for the eigenvalue c, when the solutions are extended out to the far-fietd 
potential regions. 

Inspection suggests that the solution (4.6) diverges as Y becomes large. The precise rate of diver- 
gence depends on the asymptotic behavior of the basic flow u. For Blasius flow this it is given by 
u~ 1 - (A/|Y|)e" Yl as |Y| — > while for the far wake we have u ~ 1 - (l/x 1/2 )e Y . It can be anti- 
cipated from the results of Gill and Drazin (1965), however, that the exact asymptotic behavior is not 
important and the same final result is obtained to order a m in the wave-speed at least. In order to 
illustrate the method, therefore, we choose a model asymptotic behavior at infinity. The reason for this 
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is that tiie level where the expansion breaks down can be found explicitly in terms of a (see below) 
and the analysis is greatly reduced. The Blasius flow gives a breakdown at a logarithmically large dis- 
tance away, together with the need of a logarithmically thin region there. The analysis bears a lot of 
resemblance to the hypersonic limit work of Cowley and Hall (1989) on the stability of compressible 
boundary layers. The final result is unaltered, however, and for brevity we take a model asymptotic 
behavior. The far-wake profile can be worked out explicitly also (it is a simple analogue of our sug- 
gested asymptotic behavior); assume, therefore that u - 1 - qe _r|Y| as | Y| -><*>. Then the third and 
fourth terms in the expansion (4.1a) become comparable when Y — (l/2r) ln(oc) due to the exponential 
growth of <j> 4 . A new region, II say, is introduced given by 

Yji = — r~lnoc + Y , 

2r 

where Y is of order one. Consideration of the solution in region I for large Y suggests the following 
expansions in II : 

4>ii = a I/2 <t>i + + ot 3/2 4> 3 + * * ' + a 5 ^ina ($j + 4> 2 + ...) + ■ • • , (4.6) 

These expansions are substituted into the Rayleigh equation (2.5), together with the appropriately 
expanded forms of G and T to give the following solutions (the solutions given here have already been 
matched to those in region I ) 

♦ l = A 0 (-qe' rY - c,) , (4,7a) 

$2 = ~ c 2Ao . (4.7b) 

y 

^3 = d 4 (— qe _rY — c,) J - c 3 A 0 , (4.7c) 

(qe" rY + c,) 2 

<j>, = (-qe' rY - c,) ~ . (4.7d) 

To complete the solution a potential region III say, of extent a -1 ia also required, so introduce a new 
variable Y = a -1 £ with £ order one. Consideration of (4.7a-d) as Y — > <» implies the following 
expansions in III ; 

<t>in = a I/2 <P 0 + aO, + o? n <t> 2 + • ■ • + e~ a ^(<D 0 + a. m O, + ...) + • • • . (4.8) 


The leading order solution is 
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O 0 = Be* , (4.9) 

where B is a constant. The final result is obtained by matching (4.9), for small with (4.7a-d) for 
large Y. We find 

d4 AqCj • 


which when coupled with the expression (4.6) for d 4 yields 


m 


Cl = ± 


I s ? 


dY 


(4.8) 


This result is in complete agreement with that of Gill and Drazin (1965) and it has been shown that 
long-wave disturbances are almost neutral, travel with the free-stream and have growth rates of order 
a m . In view of the dependence of T on Mach number, the wake is more stable the hotter it is, i.e. the 
higher tha Mach number. Even though the growth rates are small with respect to the free-stream, they 
become leading order effects in a frame of reference that travels with unit speed and so a classical 
weakly nonlinear approach is not likely to work (see PS2). The situation in the hypersonic limit seems 

—9 

different, however, since the growth rate decays to zero like M now. 


5. Numerical solutions. 

The previous sections presented some aspects of the analytical properties of the instability in the 
special cases of short- and long-waves. As shown the growth rates are asymptotically small and in 
order to find the maximally growing wavelengths or frequencies, as the case may be, the problem must 
be addressed numerically in a systematic search of eigenvalues. The central element of our computa- 
tions is the use of correct mean-flow profiles in the calculation of stability characteristics. Particularly, 
this becomes essential in the near-wake where the use of model, e.g. gaussian profiles, is not appropri- 
ate. The steady equations (2.3) together with (2.3BC) are solved numerically by marching forward in 
x (Smith (1974), Cebeci et al. (1979), PS1). Some results for u and T are shown in Figures 1 and 2. 
Refinement checks were made to ensure numerical accuracy. Typically 600 points are used, 300 in the 
Goldstein layer and 300 in the boundary layer. Solutions are depicted at various increasing streamwise 
locations ranging from £ = 0.2 to 2.0 (i.e. the x variation is from 0.008 to 8.0), and due to symmetry 
the region y > 0 alone is shown. At the smallest value of ^ the profile contains the Goldstein double- 
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structure quite clearly, the main part of the flow being the displaced Blasius solution. Further down- 
stream the center-line velocity increases and the profile develops into its far-wake form which is a 
gnussian perturbation from the uniform flow at infinity. The center-line temperature, on the other hand, 
decreases as the distance from the trailing edge increases with a return to the ambient temperature at 
downstream infinity. The present work concentrates on the stability of the near-wake region where the 
initial growth of disturbances takes place. It seems likely that further downstream the wake dynamics 
become nonlinear as has been observed in experimental work on incompressible wakes (Sato and 
Kuriki (I960)). 

The computational strategy is as follows. A streamwise station, is chosen and the correspond- 
ing steady-state velocity and temperature profiles are obtained. This mean flow is then used in a 
quasi-parallel manner (essentially % is a parameter) to define the compressible Rayleigh problem (2.4), 
(2.4BC). For a given Mach number (2.4) is integrated by starting at a large value of Y (usually the 
outer limit of the mean-flow calculation) where the eigenfunction is generated by its exponentially 
small behavior - this is readily found from (2.4) and the decaying solution is chosen in order to satisfy 
boundedness at infinity (the possibility of neutral oscillatory solutions as the ones found in shear 
layers, do not arise in the present problem; see later). A fourth order Rungc-Kutta integration extends 
the solution to Y = 0, where the choice of either mode I or mode II disturbances fixes a linear eigen- 
value problem which is solved by a root-finding routine. Most of the calculations address spatial stabil- 
ity whereby the frequency co = etc is taken to be real and the eigenvalues cc and c are generally com- 
plex. This stability problem corresponds to disturbances growing or decaying with downstream distance 
and seems to be more relevant that the temporal problem where disturbances grow in time at fixed 
positions (see later discussion on convective and absolute instability). Neutral eigensolutions are also 
constructed, and as has been shown by Lees and Lin (1946) these have wave-speeds equal to the value 
of the mean flow at the generalized inflection points. If these points are denoted by Y c , then they are 
given by the zeros of 


L(Y) = — 
dY 


_Uy 

T 2 


Non-singular neutral eigensolutions have wave-speeds equal to u(Y c ). In this case the Rayleigh equa- 
tion has a removable singularity at Y c and in computing neutral eigensolutions we used deformation of 


lllli Ml I II innmhi n iiiiiMiiitHiiiiii N n i n mi mindniiiiiMMi p Imnm ihmi'i i ,«li 



- 19 - 


the integration into the complex plane below the singularity. By numerical evaluation of L(Y) it can 
be concluded that the wake-flows under consideration have two generalized inflection points which are 
symmetrically placed about the wake center-line. This situation persists at arbitrarily large Mach 
numbers also, even though Y c now gets pushed towards the outer edges of the wake. In Figure 3 we 
show a collection of L(Y) at different increasing Mach numbers and for the half-space Y > 0; the 
profile used is at £ = 0.8 but the behavior is qualitatively the same at other wake-stations also. This in 
turn means that, for a given x-station and Mach number, all neutral waves have the same wave-speed. 


6. Results and discussion. 

Before the presentation of stability results it is useful to consider the physical nature of the possi- 
ble modes. A lot of information can be found by consideration of the eigenfunctions at infinity. Our 
discussion follows the classification of instability waves in shear layers by Jackson and Grosch (1989). 
The condition imposed on the disturbance is that it remains bounded as Y -) «, For large Y, there- 
fore, (2.4) gives the following behavior for 2D disturbances (3D ones form a simple extension, and we 
only want to illustrate things here) : 

p~e ±nAY , A 2 = 1 — (l-c) 2 M z . (6.1) 

The zeros of of A define the sonic lines where the disturbance at infinity neither grows nor decays. 
These are given by 

C±=,± ^ • 

It can be seen from this equation that the two sonic lines do not cross for any Mach number, and 
further more c + represents a wave with speed higher than the free-stream value and is therefore pre- 
cluded on physical grounds. If A 2 is positive, the negative sign in the exponential in (6.1) is taken and 
the disturbance vanishes at infinity - the perturbation is termed subsonic since it has characteristics 
similar to those of incompressible stability. If A 2 is negative, however, the solutions are oscillatory at 
infinity where they take the form of outgoing or incoming waves. Such neutral waves are singular in 
that they do not satisfy the Lees and Lin regularity condition, but they are physically relevant as zero 
viscosity limit solutions of the "Orr-Sommerfeld" problem (Lees and Lin (1946), Jackson and Grosch 
(1989)). In the present wake problem the possibility of A 2 being negative does not arise, and to see 
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this wc construct the neutral wave-speed from the Lees and Lin generalized inflection point criterion 
and plot it together with the equation of the sonic line c_. A typical situation is depicted in Figure 4 for 
two wake stations ^ = 0.1 and ^ = 1.2. The conclusion is, therefore, that subsonic boundary conditions 
at infinity are adequate. 

As mentioned earlier, we chose to carry out a spatial stability analysis which is, in general, more 
pertinent to wake-flows than a temporal one. Before presenting results it is worth mentioning that the 
accuracy of computed eigenvalues was checked by doubling the number of points in the mean-flow 
computation. The finding is that the representation of the mean-flow by 600 points produces more than 
graphical accuracy in the results. 

First we address the convective/absolute instability question. Tn Figure 5a we show a collection 
of mode I spatial stability curves depicting the variation of -oq (the spatial growth-rate) with frequency 
and for the basic flow profiles at the streamwise locations 
£ = 0.8 , 0.6 , 0.4 , 0.2 , 0.08 , 0.06 , 0.04 , 0X)2 , 0.01. The Mach number is 0.1 and the correspond- 
ing wave-number curves ( a,, vs CO ) are given in Figure 5b. Figure 5a shows that as we get closer to 
the trailing edge the maximum growth rate increases sharply and the frequency band of unstable waves 
decreases. This behavior is very similar to that discussed in PS1 for incompressible near-wakes. The 
fact that the instability is changing from convective to absolute, can be seen from Figure 5b. The 
curves become steeper as £ decreases and by the stability at the final wake-station ^ = 0.01 is charac- 
terized by an infinite gradient. The value of d(o/da. r is tending to zero, which is the condition for abso- 
lute instability (see Huerre and Monkewitz (1985)). It is found, therefore, that there exists a region of 
absolute instability in the near wake for small Mach numbers; it is worth noting that the extent of this 
region is very small, less than 0.00000 lxL in dimensional terms. 

The effect of Mach number on the presence of absolute instability regions is investigated next. 
The same wake-stations as in Figures 5 are used but the Mach number is now 3.0. Figures 6a and 6b 
show the variation of spatial growth-rates and wave-numbers with frequency, for the various wake- 
stations in question. First, note that in absolute terms the maximum growth rates decrease by more than 
a factor of ten and the unstable frequency band decreases by a factor of two about. In general terms, 
therefore, the effect of a higher Mach number is highly stabilizing (see results later also for more quan- 
titative conclusions) The sharp features and singularity formation of Figures 5a, b are no longer present. 
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In fact, as Figure 6a shows, the spatial growth rate curves converge to a limiting one as £ is decreased 
(the change between the curves at £ = 0.02 and £ = 0.01 is almost negligible). The fact that no singu- 
larity forms is clarified further by Figure 6b; the variation of a, with frequency tends to a limiting 
curve with a finite slope. 

The previous results addressed the effect of the basic-flow development on its stability charac- 
teristics. In what follows we pick two typical wake-stations and carry out a more detailed stability cal- 
culation. The two profiles used correspond to the values £ = 0.1 and £ = 1.2. They differ in that the 
former is a near-wake profile ( x = 0.001 ) while the latter is a far-wake one ( x = 1.728). Their 
respective center-line velocities are 0.0772 and 0.72. Physically the near-wake still retains many of the 
characteristics of the upstream boundary-layer flow while in the far-wake the double structures seen 
upstream are no longer present. 

Figures 7a and 7b depict the effect of Mach number on the wake-stability at £ = 0.1. Mode I 
waves are computed at Mach numbers 0, 0.5, 1, 1.5, 2, 2.5 and 3. A marked increase in stability is 
observed with the maximum growth rate dropping from approximately 0.27 at zero Mach number to 
0.03 at a Mach number of 3. The frequency range of unstable waves also decreases. The value of the 
neutral wave-number ( a ) also decreases with Mach number as can be seen from Figure 7b; the neu- 
tral a is the value at the ends of the depicted curves since these points correspond to the almost 
attained neutral points on the spatial stability curves 7a. The curves go from right to left as the Mach 
number increases. This calculation was carried further to the higher Mach numbers of 4, 5, 6, 7, 8 and 
9. The results are shown in Figures 8a and 8b. The continued increase of stability is clearly demon- 
strated and at a Mach number of 9 the maximum growth rate is tiny, about 0.002. Figure 8b depicts 
the variation of 0^ with O) for the two Mach numbers of 4 and 9 only; the curves for intermediate 
Mach numbers are almost coincident but their end-points are at lower values of a r . The curve for 
M = 4 is deliberately not continued to the origin since we are mainly interested at the endpoints. The 
conclusion that can be drawn from this calculation is that the neutral a is decreasing to zero with 
Mach number. The situation is quite different at a far-wake station as is shown later. 

Figures 9a and 9b describe the stability characteristics of mode II waves at the station £ = 0.1 
and at Mach numbers of 0, 0.5 and 1.0. Comparison with Figures 7a, b shows that mode II is much 
more stable than mode I. For example at M = 1 the maximum growth rate of mode I waves is 
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approximately 0.15 while that for a mode II wave at the same Mach number is only 0.0025; the 

corresponding ranges of unstable frequencies are 0 to 0.32 and 0.006 respectively. Further more, at the 

wake-station under consideration, the maximum growth rate of the mode II waves increase slightly for 
oblique waves, but for a given Mach number they are still much lower than the corresponding 2D 
mode l growth rates. Figure 9b shows that the neutral wave-number is again decreasing with Mach 
number. 

In Figures 10a and 10b we show the effect of obliqueness on the instability. In particular, mode 
I waves are computed at a Mach number of 3 and inclined at angles of 0, 20, 30, 45, 60 and 75 

degrees respectively to the direction of the basic flow. The maximum growth rates are seen to increase 

as the angle increases, but the neutral frequencies decrease with angle. The final computed angle of 75 
degrees has the highest maximum growth rate which is almost twice that of the 2D wave. When the 
angle increases further the curve begins to come down. For small wave-propagation angles the fre- 
quency where the maximum growth rate is achieved increases by a small amount from the 2D one, but 
when the angle increases from 60 to 75 degrees the maximum growth rate frequency shifts to a lower 
value. The corresponding variation of a r with 0) is shown in Figure' 10b. The curves move from right 
to left as the angle increases. 

In Figures 11a and I lb we present stability characteristics of a typical far-wake station that has 
^ = 1.2. Mode I waves are computed and growth rate and wave-number curves are shown for the 
Mach numbers 0, 1,2, 3, 4, 5, 6 and 7. The qualitative nature of the results is very similar to those of 
the near-wake station ^ = 0.1. The maximum growth rates are much lower, however (e.g. at M = 1 the 
maximum growth rates are approximately 0.018 and 0.16 respectively). The main quantitative 
difference between the two stations is that the neutral value of a for the far-wake profile is approach- 
ing a non-zero constant as the Mach number increases, whereas the near-wake value was decreasing to 
zero. This apparent difference in the physics is not too surprising in view of the completely different 
structure of the basic flow in the near- and far-wake respectively. 

7. Conclusions. 

The present study is concerned with the understanding of the physical instability mechanisms that 
are at play in free shear flows and in particular symmetric wakes. The main emphasis was given to the 
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instability of near-wakes. It was shown that these mechanisms can only be captured correctly if the 
correct basic flow is used. In fact our results suggest that the hypersonic limit of far-wakes is very 
different from that of near-wakes (this problem is currently being analyzed by the author), and the use 
of a model profile (e.g. gaussian distribution) in an attempt to capture the near-wake dynamics is 
incorrect. Using our computed basic-flow we showed that for low enough Mach numbers there exists 
a small region of absolute instability very near the trailing edge which disappears as the Mach number 
increases. The maximum spatial growth rates available increase as the trailing edge is approached, 
irrespective of the Mach number. Further more mode I waves are the most dangerous in the near-wake 
and the instability (maximum growth-rate) can be almost doubled if the wave is three-dimensional. If 
a temporal theory is carried out instead, some different behavior results. For example, as Chen et al. 
(1989) report, a far-wake profile modeled by the gaussian distribution u = 1 - 0.692e 069315Y together 
with the corresponding temperature variation, gives the highest temporal growth rate for mode I waves 
which are two dimensional, and at all Mach numbers. In Figure 12 we show the spatial stability of this 
profile, which shows that the growth rate is larger for 3D waves inclined at an angle of about 60 
degrees. Such discrepancies between spatial and temporal calculations can, therefore, be very mislead- 
ing in practical applications. 
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Figure 4 Mode classification diagram. Solid curve is the sonic line, the other curves show the variation 
of the neutral wave-speed with Mach number for the basic (low at ^ = 0.1 and £ = 1.2. 
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M = 0. 1 £ = .8,.6,.4,.2,.1 ,.08,.06,.04,.02,.01 



Figure 5a Spatial stability of proper compressible wake profiles at various streamwise locations. Mach 
number M=0.1. The graphs show the variation of the spatial growth rate -otj with frequency co. The 
curves correspond to the streamwise locations £ = 0.8, 0.6, 0.4, 0.2, 0.1, 0.08, 0.06, 0.04, 0.02 and 


0j 
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Figure 6a Spatial stability of proper compressible wake profiles at various streamwise locations. Mach 
number M=3.0. The graphs show the variation of the spatial growth rate -a, with frequency to. The 
curves correspond to the streamwise locations £ = 0.8, 0.6, 0.4, 0.2, 0.1, 0.08, 0.06, 0.04, 0.02 and 
0j 
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Figure 6b Same as Figure 6a but the variation of a T with o> is now shown. 
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Figure 10a Near wake profile £ = 0.1. 3D mode I waves. S 



Mach number of 3.0. The angles of propagation are 0 = 0, 20 
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Figure 10b Same as Figure 10a but variation of a, with o. 
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Figure 11a Far wake profile = 1.2. Spatial stability growth rate curves of -alphnj vs w for M=0, I, 
2. 3, 4, 5, 6 and 7. 
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Figure lib Same as 11a but variation of a,, with w. 
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Figure 12 Gaussian model profile - 3D spatial stability at angles 0= 0, 30, 60, 75. 
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